Structural Impact Assessment of Cytochrome P450 2A13 Polymorphisms Using Molecular Dynamics Simulations.

One of the members of CYP, a monooxygenase, CYP2A13 is involved in the metabolism of nicotine, coumarin, and tobacco-specific nitrosamine. Genetic polymorphisms have been identified in CYP2A13, with reported loss or reduction in enzymatic activity in CYP2A13 allelic variants. This study aimed to unravel the mechanism underlying the diminished enzymatic activity of CYP2A13 variants by investigating their three-dimensional structures through molecular dynamics (MD) simulations. For each variant, MD simulations of 1000 ns were performed, and the obtained results were compared with those of the wild type. The findings indicated alterations in the interaction with heme in CYP2A13.4, .6, .8, and .9. In the case of CYP2A13.5, observable effects on the helix structure related to the interaction with the redox partner were identified. These conformational changes were sufficient to cause a decrease in enzyme activity in the variants. Our findings provide valuable insights into the molecular mechanisms associated with the diminished activity in the CYP2A13 polymorphisms.

CYP2A6 associates with respiratory disease risk and younger age of diagnosis: a phenome-wide association Mendelian Randomization study.

CYP2A6, a genetically variable enzyme, inactivates nicotine, activates carcinogens, and metabolizes many pharmaceuticals. Variation in CYP2A6 influences smoking behaviors and tobacco-related disease risk. This phenome-wide association study examined associations between a reconstructed version of our weighted genetic risk score (wGRS) for CYP2A6 activity with diseases in the UK Biobank (N = 395 887). Causal effects of phenotypic CYP2A6 activity (measured as the nicotine metabolite ratio: 3'-hydroxycotinine/cotinine) on the phenome-wide significant (PWS) signals were then estimated in two-sample Mendelian Randomization using the wGRS as the instrument. Time-to-diagnosis age was compared between faster versus slower CYP2A6 metabolizers for the PWS signals in survival analyses. In the total sample, six PWS signals were identified: two lung cancers and four obstructive respiratory diseases PheCodes, where faster CYP2A6 activity was associated with greater disease risk (Ps < 1 × 10-6). A significant CYP2A6-by-smoking status interaction was found (Psinteraction < 0.05); in current smokers, the same six PWS signals were found as identified in the total group, whereas no PWS signals were found in former or never smokers. In the total sample and current smokers, CYP2A6 activity causal estimates on the six PWS signals were significant in Mendelian Randomization (Ps < 5 × 10-5). Additionally, faster CYP2A6 metabolizer status was associated with younger age of disease diagnosis for the six PWS signals (Ps < 5 × 10-4, in current smokers). These findings support a role for faster CYP2A6 activity as a causal risk factor for lung cancers and obstructive respiratory diseases among current smokers, and a younger onset of these diseases. This research utilized the UK Biobank Resource.

Characterization of CYP2B6 and CYP2A6 Pharmacogenetic Variation in Sub-Saharan African Populations.

Genetic variation in CYP2B6 and CYP2A6 is known to impact interindividual response to antiretrovirals, nicotine, and bupropion, among other drugs. However, the full catalogue of clinically relevant pharmacogenetic variants in these genes is yet to be established, especially across African populations. This study therefore aimed to characterize the star allele (haplotype) distribution in CYP2B6 and CYP2A6 across diverse and understudied sub-Saharan African (SSA) populations. We called star alleles from 961 high-depth full genomes using StellarPGx, Aldy, and PyPGx. In addition, we performed CYP2B6 and CYP2A6 star allele frequency comparisons between SSA and other global biogeographical groups represented in the new 1000 Genomes Project high-coverage dataset (n = 2,000). This study presents frequency information for star alleles in CYP2B6 (e.g., *6 and *18; frequency of 21-47% and 2-19%, respectively) and CYP2A6 (e.g., *4, *9, and *17; frequency of 0-6%, 3-10%, and 6-20%, respectively), and predicted phenotypes (for CYP2B6), across various African populations. In addition, 50 potentially novel African-ancestry star alleles were computationally predicted by StellarPGx in CYP2B6 and CYP2A6 combined. For each of these genes, over 4% of the study participants had predicted novel star alleles. Three novel star alleles in CYP2A6 (*54, *55, and *56) and CYP2B6 apiece, and several suballeles were further validated via targeted Single-Molecule Real-Time resequencing. Our findings are important for informing the design of comprehensive pharmacogenetic testing platforms, and are highly relevant for personalized medicine strategies, especially relating to antiretroviral medication and smoking cessation treatment in Africa and the African diaspora. More broadly, this study highlights the importance of sampling diverse African ethnolinguistic groups for accurate characterization of the pharmacogene variation landscape across the continent.

Associating CYP2A6 structural variants with ovarian and lung cancer risk in the UK Biobank: replication and extension.

CYP2A6 is a polymorphic enzyme that inactivates nicotine; structural variants (SVs) include gene deletions and hybrids with the neighboring pseudogene CYP2A7. Two studies found that CYP2A7 deletions were associated with ovarian cancer risk. Using their methodology, we aimed to characterize CYP2A6 SVs (which may be misidentified by prediction software as CYP2A7 SVs), then assess CYP2A6 SV-associated risk for ovarian cancer, and extend analyses to lung cancer. An updated reference panel was created to impute CYP2A6 SVs from UK Biobank array data. Logistic regression models analyzed the association between CYP2A6 SVs and cancer risk, adjusting for covariates. Software-predicted CYP2A7 deletions were concordant with known CYP2A6 SVs. Deleterious CYP2A6 SVs were not associated with ovarian cancer (OR = 1.06; 95% CI: 0.80-1.37; p = 0.7) but did reduce the risk of lung cancer (OR = 0.44; 95% CI: 0.29-0.64; p < 0.0001), and a lung cancer subtype. Replication of known lung cancer associations indicates the validity of array-based SV analyses.

CYP2A6 Activity and Cigarette Consumption Interact in Smoking-Related Lung Cancer Susceptibility.

Cigarette smoke, containing both nicotine and carcinogens, causes lung cancer. However, not all smokers develop lung cancer, highlighting the importance of the interaction between host susceptibility and environmental exposure in tumorigenesis. Here, we aimed to delineate the interaction between metabolizing ability of tobacco carcinogens and smoking intensity in mediating genetic susceptibility to smoking-related lung tumorigenesis. Single-variant and gene-based associations of 43 tobacco carcinogen-metabolizing genes with lung cancer were analyzed using summary statistics and individual-level genetic data, followed by causal inference of Mendelian randomization, mediation analysis, and structural equation modeling. Cigarette smoke-exposed cell models were used to detect gene expression patterns in relation to specific alleles. Data from the International Lung Cancer Consortium (29,266 cases and 56,450 controls) and UK Biobank (2,155 cases and 376,329 controls) indicated that the genetic variant rs56113850 C>T located in intron 4 of CYP2A6 was significantly associated with decreased lung cancer risk among smokers (OR = 0.88, 95% confidence interval = 0.85-0.91, P = 2.18 × 10-16), which might interact (Pinteraction = 0.028) with and partially be mediated (ORindirect = 0.987) by smoking status. Smoking intensity accounted for 82.3% of the effect of CYP2A6 activity on lung cancer risk but entirely mediated the genetic effect of rs56113850. Mechanistically, the rs56113850 T allele rescued the downregulation of CYP2A6 caused by cigarette smoke exposure, potentially through preferential recruitment of transcription factor helicase-like transcription factor. Together, this study provides additional insights into the interplay between host susceptibility and carcinogen exposure in smoking-related lung tumorigenesis. SIGNIFICANCE: The causal pathway connecting CYP2A6 genetic variability and activity, cigarette consumption, and lung cancer susceptibility in smokers highlights the need for behavior modification interventions based on host susceptibility for cancer prevention.

Cigarette Smoking and Intracranial Aneurysms: A Pilot Analysis of SNPs in the CYP2A6 Gene in the Italian Population.

BACKGROUND: Cigarette smoking is a modifiable risk factor associated with formation and rupture of intracranial aneurysms (IAs). Cytochrome P450 2A6 (CYP2A6) is the main enzyme implied in catabolism of nicotine and xenobiotics, giving rise to oxidative stress products. Our study investigated the associations between specific single-nucleotide polymorphisms (SNPs) in the CYP2A6 gene and the presence of sporadic IAs in a cluster of Italian patients, as well as their rupture regarding cigarette smoking habit. METHODS: Three hundred and thirty-one Italian patients with sporadic IAs were recruited in a single institution. We recorded data on clinical onset with subarachnoid hemorrhage (SAH) and smoking habit. Genetic analysis was performed with a standard procedure on peripheral blood samples: CYP2A6 ∗1B2, CYP2A6 ∗2, and CYP2A6 ∗14 SNPs were analyzed in the study group along with 150 healthy control subjects. Statistical analysis was conducted according to genetic association study guidelines. RESULTS: In the patient cohort, the frequency of aSAH was significantly higher in current smokers (P < 0.001; OR=17.45), regardless of the pattern of CYP2A6 SNPs. There was a correlation between IA rupture and cigarette smoking in patients with the heterozygous CYP2A6 ∗1B2 allele (P < 0.001; OR=15.47). All patients carrying the heterozygous CYP2A6 ∗14 allele had an aSAH event (100%), regardless of smoking habit, although this correlation was not statistically significant (P = 1). CONCLUSIONS: According to our findings, a cigarette smoker carrying a fully active CYP2A6 enzyme (heterozygous ∗1B2 allele) may have an increased risk of IA rupture compared to those with functionally less active variants: further investigation on a larger sample is needed to verify this result. The role of the heterozygous CYP2A6 ∗14 allele in aSAH is yet to be clarified.

Recurrent CYP2A6 gene mutation in biphasic hyalinizing psammomatous renal cell carcinoma: Additional support of three cases.

Biphasic hyalinizing psammomatous renal cell carcinoma (BHP RCC) with NF2 gene mutations is a newly described provisional category of renal cell carcinoma (RCC). Here we described three additional cases of BHP RCC with CYP2A6 gene mutation besides NF2 gene. The carcinomas were predominantly unencapsulated, and two of them had a rounded, nodular interface with the native kidney while one had perirenal adipose tissue invasion. Histopathologically, all neoplasms had a characteristic biphasic appearance of smaller cells clustering around basement membrane material within larger acini, forming pseudorosettes or a glomeruloid pattern. The smaller cells were focally spindle-shaped in two carcinomas. Psammoma bodies were shown in two carcinomas. Cellular necrosis and perineural invasion was identified in one case. Immunohistochemically, Vimentin, EMA, P504s were extensively expressed while RCC and CD10 were only expressed in larger cells. CK7 was positive in one tumor. CYP2A6 gene mutation (CYP2A6 NM_000762.6: exon4:c.A580G:p.K194E) was revealed in three tumors by Whole-genome exome sequencing, which was further confirmed by Sanger sequencing. Only one case harbored a somatic termination mutation in NF2 gene. NF2 promoter methylation was observed in the other two cases. Clinically, one patient died of disease with widespread bone metastases confirmed by biopsy at the ninth month after surgery but the other two patients had no evidence of recurrence or metastases (follow-up period 9-90 months). Our findings validated previously described clinicopathological features and NF2 gene mutation or promoter methylation of BHP RCC. In addition, we reported different IHC pattern of BHP RCC and further revealed the recurrent CYP2A6 genetic alteration.

Genotyping, characterization, and imputation of known and novel CYP2A6 structural variants using SNP array data.

CYP2A6 metabolically inactivates nicotine. Faster CYP2A6 activity is associated with heavier smoking and higher lung cancer risk. The CYP2A6 gene is polymorphic, including functional structural variants (SV) such as gene deletions (CYP2A6*4), duplications (CYP2A6*1 × 2), and hybrids with the CYP2A7 pseudogene (CYP2A6*12, CYP2A6*34). SVs are challenging to genotype due to their complex genetic architecture. Our aims were to develop a reliable protocol for SV genotyping, functionally phenotype known and novel SVs, and investigate the feasibility of CYP2A6 SV imputation from SNP array data in two ancestry populations. European- (EUR; n = 935) and African- (AFR; n = 964) ancestry individuals from smoking cessation trials were genotyped for SNPs using an Illumina array and for CYP2A6 SVs using Taqman copy number (CN) assays. SV-specific PCR amplification and Sanger sequencing was used to characterize a novel SV. Individuals with SVs were phenotyped using the nicotine metabolite ratio, a biomarker of CYP2A6 activity. SV diplotype and SNP array data were integrated and phased to generate ancestry-specific SV reference panels. Leave-one-out cross-validation was used to investigate the feasibility of CYP2A6 SV imputation. A minimal protocol requiring three Taqman CN assays for CYP2A6 SV genotyping was developed and known SV associations with activity were replicated. The first domain swap CYP2A6-CYP2A7 hybrid SV, CYP2A6*53, was identified, sequenced, and associated with lower CYP2A6 activity. In both EURs and AFRs, most SV alleles were identified using imputation (>70% and >60%, respectively); importantly, false positive rates were <1%. These results confirm that CYP2A6 SV imputation can identify most SV alleles, including a novel SV.

Influence of CYP2A6 Genetic Variation, Nicotine Dependence Severity, and Treatment on Smoking Cessation Success.

INTRODUCTION: Genetic variation in Cytochrome P450 2A6 (CYP2A6), the major nicotine metabolizing enzyme, is associated with nicotine dependence and smoking cessation. Nicotine dependence severity also predicts smoking cessation. Our goals were to determine how CYP2A6 variation and nicotine dependence alter smoking cessation, and whether dependence could refine CYP2A6-based treatment recommendations. AIMS AND METHODS: Adult smokers treated for 12 weeks with placebo, nicotine patch, or varenicline (NCT01314001) were grouped as CYP2A6 normal (n = 567) or slow (n = 432) nicotine metabolizers based on a CYP2A6 weighted genetic risk score. Fagerström test for nicotine dependence scores were measured at baseline and biochemically verified smoking cessation was assessed at end of treatment. RESULTS: Dependence neither mediated nor moderated an association between CYP2A6 variation and smoking cessation overall, within any treatment arm, or after stratifying by ancestry (n = 591 European, n = 408 African ancestry) or sex (n = 444 women, n = 555 men). In within-treatment analyses, the mediation effect odds ratio (OR) ranged from 0.95 to 1.00 and the bias-corrected 95% confidence interval contained 1. Moderation (i.e. interaction) effect ORs ranged from 0.88 to 1.61 (p = .397-.828). For CYP2A6 normal metabolizers, quit rates on varenicline were similar for those with high (41.1%) and low (43.4%) dependence, while quit rates were lower for those with high versus low dependence on both patch (16.5 vs. 29.7%) and placebo (8.9 vs. 18.5%). CYP2A6 slow metabolizers with high versus low dependence had lower quit rates in all three treatment arms. CONCLUSIONS: Although nicotine dependence severity neither mediated nor moderated CYP2A6 associations with smoking cessation, incorporating information on dependence may optimize the choice of smoking cessation treatment aid in CYP2A6 normal and slow metabolizers. IMPLICATIONS: Variation in CYP2A6 and nicotine dependence severity alter smoking cessation success. Our findings suggest that while nicotine dependence severity is unlikely to mediate or moderate CYP2A6 associations with cessation, incorporating patient information on both CYP2A6 and nicotine dependence severity may lead to improved smoking cessation strategies.

Newly Identified Tree Shrew Cytochrome P450 2A13 is Expressed in Liver and Lung and Encodes a Functional Drug-Metabolizing Enzyme Similar to Dog Cytochrome P450 2A13 and Pig Cytochrome P450 2A19.

The tree shrew, a non-rodent primate-like species, is used in various fields of biomedical research, including hepatitis virus infection, myopia, depression, and toxicology. Recent genome analysis found that the numbers of cytochrome P450 (P450 or CYP) genes are similar in tree shrews and humans and their sequence identities are high. Although the P450s are a family of important drug-metabolizing enzymes, they have not yet been fully investigated in tree shrews. In the current study, tree shrew CYP2A13 cDNA was isolated from liver, and its characteristics were compared with those of pig, dog, and human CYP2As. Tree shrew CYP2A13 amino acid sequences were highly identical (87-92%) to the human CYP2As and contained sequence motifs characteristic of P450s. Phylogenetic analysis revealed that tree shrew CYP2A13 was more closely related to human CYP2As than to rat CYP2As, similar to dog and pig CYP2As. Among the tissue types analyzed, tree shrew CYP2A13 mRNA was preferentially expressed in liver and lung, similar to dog CYP2A13 mRNA, whereas dog CYP2A25 and pig CYP2A19 mRNAs were predominantly expressed in liver. Tree shrew liver microsomes and tree shrew CYP2A13 proteins heterologously expressed in Escherichia coli catalyzed coumarin 7-hydroxylation and phenacetin O-deethylation, just as human, dog, and pig CYP2A proteins and liver microsomes do. These results demonstrate that tree shrew CYP2A13 is expressed in liver and lung and encodes a functional drug-metabolizing enzyme. SIGNIFICANCE STATEMENT: Novel tree shrew cytochrome P450 2A13 (CYP2A13) was identified and characterized in comparison with human, dog, and pig CYP2As. Tree shrew CYP2A13 isolated from liver had high sequence identities and close phylogenetic relationships to its human homologs and was abundantly expressed in liver and lung at the mRNA level. Tree shrew CYP2A13 metabolized coumarin and phenacetin, human selective CYP2A6 and CYP2A13 substrates, respectively, similar to dog and pig CYP2As, and is a functional drug-metabolizing enzyme likely responsible for drug clearances.

Mining of molecular insights of CYP2A6 and its variants complex with coumarin (CYP2A6*-coumarin) using molecular dynamics simulation.

CYP2A6 is a very important enzyme that plays a crucial role in nicotine compounds and is responsible for the metabolism of more than 3% drugs of total metabolized drugs by the CYP family and reported as one of very important pharmacogenes. CYP2A6 is highly polymorphic in nature and reported with more than 40 variants, most of these variants are SNPs originated and population specific. It has been well observed and reported that the presence of these population-specific non-synonymous SNPs in CYP2A6 alters the rate of drug metabolism and as a functional consequence, drugs produce an abnormal response. Though genomics and pharmacogenomics studies are there, very less is known about the structural effects of these SNPs on molecular-interaction and folding of CYP2A6. To fill the knowledge gap, SNPs based four variants, i.e., CYP2A6*2, CYP2A6*18, CYP2A6*21, and CYP2A6*35, which are frequently reported in the South Asian population, were considered for the study. Coumarin (DB04665), a well reported drug, is considered as a model substance, and the effect of all four variants on 'CYP2A6*-coumarin' complex was studied. MD simulation-based analysis (at 200 ns) was performed and comparative analysis with respect to wild type 'CYP2A6-coumarin' complex was done. Though observation didn't find any global effect on complete complex but found some crucial minor-local alteration in interaction and folding process. It is assumed that the change due to SNPs in the single amino acid did not bring global change in physiochemical properties of CYP2A6* but caused local-trivial changes which are very crucial for its metabolic activity.Communicated by Ramaswamy H. Sarma.

Impact of Genetic Variants in the Nicotine Metabolism Pathway on Nicotine Metabolite Levels in Smokers.

BACKGROUND: Nicotine metabolism is a major factor in nicotine dependence, with approximately 70% to 80% of nicotine metabolized to cotinine in Caucasians. Cotinine formation is catalyzed primarily by CYP2A6, which also converts cotinine to trans-3'-hydroxycotinine (3HC). The goal of the present study was to examine the effects of CYP2A6 deficiency on nicotine metabolism profiles in vivo and the importance of genetic variants in nicotine-metabolizing enzyme genes on urinary nicotine metabolites levels. METHODS: Urine samples from 722 smokers who participated in the Singapore Chinese Health Study were analyzed using UPLC-MS/MS to detect nicotine and eight of its urinary metabolites, and a total of 58 variants in 12 genes involved in nicotine metabolism were investigated in 475 of these subjects with informative genotyping data. RESULTS: Urine samples stratified by the ratio of 3HC/cotinine exhibited a 7-fold increase in nicotine-N'-oxide, a 6-fold increase in nicotine-Glucuronide (Gluc), and a 5-fold decrease in 3HC-Gluc when comparing the lower versus upper 3HC/cotinine ventiles. Significant (P < 0.0001) associations were observed between functional metabolizing enzyme genotypes and levels of various urinary nicotine metabolites, including CYP2A6 genotype and levels of nicotine, nicotine-Gluc, nicotine-N'-oxide and 3HC, UGT2B10 genotype and levels of cotinine, nicotine-Gluc and cotinine-Gluc, UGT2B17 genotype and levels of 3HC-Gluc, FMO3 genotype and levels of nicotine-N'-oxide, and CYP2B6 genotype and levels of nicotine-N'-oxide and 4-hydroxy-4-(3-pyridyl)-butanoic acid. CONCLUSIONS: These data suggest that several pathways are important in nicotine metabolism. IMPACT: Genotype differences in several nicotine-metabolizing enzyme pathways may potentially lead to differences in nicotine dependence and smoking behavior and cessation.

Accuracy and applications of sequencing and genotyping approaches for CYP2A6 and homologous genes.

OBJECTIVES: We evaluated multiple genotyping/sequencing approaches in a homologous region of chromosome 19, and investigated associations of two common 3'-UTR CYP2A6 variants with activity in vivo. METHODS: Individuals (n = 1704) of European and African ancestry were phenotyped for the nicotine metabolite ratio (NMR), an index of CYP2A6 activity, and genotyped/sequenced using deep amplicon exon sequencing, SNP array, genotype imputation and targeted capture sequencing. Amplicon exon sequencing was the gold standard to which other methods were compared within-individual for CYP2A6, CYP2A7, CYP2A13, and CYP2B6 exons to identify highly discordant positions. Linear regression models evaluated the association of CYP2A6*1B and rs8192733 genotypes (coded additively) with logNMR. RESULTS: All approaches were ≤2.6% discordant with the gold standard; discordant calls were concentrated at few positions. Fifteen positions were discordant in >10% of individuals, with 12 appearing in regions of high identity between homologous genes (e.g. CYP2A6 and CYP2A7). For six, allele frequencies in our study and online databases were discrepant, suggesting errors in online sources. In the European-ancestry group (n = 935), CYP2A6*1B and rs8192733 were associated with logNMR (P < 0.001). A combined model found main effects of both variants on increasing logNMR. Similar trends were found in those of African ancestry (n = 506). CONCLUSION: Multiple genotyping/sequencing approaches used in this chromosome 19 region contain genotyping/sequencing errors, as do online databases. Gene-specific primers and SNP array probes must consider gene homology; short-read sequencing of related genes in a single reaction should be avoided. Using improved sequencing approaches, we characterized two gain-of-function 3'-UTR variants, including the relatively understudied rs8192733.

Nicotine Metabolism Predicted by CYP2A6 Genotypes in Relation to Smoking Cessation: A Systematic Review.

INTRODUCTION: Identifying genetic factors associated with smoking cessation could inform precision cessation interventions. Of major interest is genetic variation in nicotine metabolism, largely predicted by CYP2A6 variations. AIMS AND METHODS: We conducted a systematic literature review to summarize the population-based evidence of the association between CYP2A6 and smoking cessation. In the 12 studies meeting the inclusion criteria, the known functional metabolic effect of CYP2A6 variants was used to classify nicotine metabolism as normal (>75% metabolic activity), intermediate (50.1%-75% activity), slow (25%-50% activity), and poor (<25% activity). Summary odds ratios of smoking cessation were calculated across metabolic groups, stratified by ancestry and whether participants received pharmacotherapy or placebo/no treatment. RESULTS: Among untreated people of European ancestry (n = 4 studies), those with CYP2A6 reduced metabolism were more likely to quit smoking than those with normal metabolism (Summary OR = 2.05, 95% CI 1.23 to 3.42) and the likelihood of cessation increased as nicotine metabolism decreased. Nicotine replacement therapy attenuated the association at end-of-treatment, while bupropion modified the association such that intermediate/slow metabolizers were less likely to quit than normal metabolizers (Summary OR = 0.86, 95% CI 0.79 to 0.94). Among untreated Asian people (n = 3 studies), results differed compared with those with European ancestry: those with slow metabolism were less likely to have quit smoking than normal metabolizers (Summary OR = 0.52, 95% CI 0.38 to 0.71). Evidence for people of African ancestry (n = 1 study) suggested the CYP2A6 association with cessation may differ compared with those of European ancestry. CONCLUSIONS AND IMPLICATIONS: Most studies included in this review were of European ancestry populations; these showed slower nicotine metabolism was associated with increased likelihood of smoking cessation in a dose-related manner. Pharmacotherapy appeared to attenuate or modify this association among people of European ancestry, but it is unclear whether the change in the association remains consistent after treatment ceases. This finding has implications for precision medicine cessation interventions. Based on only a few studies of people of Asian or African ancestry, the association between CYP2A6 variants and cessation may differ from that observed among those of European ancestry, but more evidence is needed.

Functional characterization of novel rare CYP2A6 variants and potential implications for clinical outcomes.

CYP2A6 activity, phenotyped by the nicotine metabolite ratio (NMR), is a predictor of several smoking behaviors, including cessation and smoking-related disease risk. The heritability of the NMR is 60-80%, yet weighted genetic risk scores (wGRSs) based on common variants explain only 30-35%. Rare variants (minor allele frequency <1%) are hypothesized to explain some of this missing heritability. We present two targeted sequencing studies where rare protein-coding variants are functionally characterized in vivo, in silico, and in vitro to examine this hypothesis. In a smoking cessation trial, 1687 individuals were sequenced; characterization measures included the in vivo NMR, in vitro protein expression, and metabolic activity measured from recombinant proteins. In a human liver bank, 312 human liver samples were sequenced; measures included RNA expression, protein expression, and metabolic activity from extracted liver tissue. In total, 38 of 47 rare coding variants identified were novel; characterizations ranged from gain-of-function to loss-of-function. On a population level, the portion of NMR variation explained by the rare coding variants was small (~1%). However, upon incorporation, the accuracy of the wGRS was improved for individuals with rare protein-coding variants (i.e., the residuals were reduced), and approximately one-third of these individuals (12/39) were re-assigned from normal to slow metabolizer status. Rare coding variants can alter an individual's CYP2A6 activity; their integration into wGRSs through precise functional characterization is necessary to accurately assess clinical outcomes and achieve precision medicine for all. Investigation into noncoding variants is warranted to further explain the missing heritability in the NMR.

Polymorphism in cytochrome P4502A6 reduces the risk to head and neck cancer and modifies the treatment outcome.

The present case-control study consisting of 1300 cases of head and neck squamous cell carcinoma (HNSCC) and the equal number of controls aimed to investigate the association of functionally important polymorphisms in cytochrome P4502A6 (CYP2A6*1B, CYP2A6*4C, CYP2A6*9-rs28399433) with HNSCC and the treatment response in cases receiving a combination of chemotherapy/radiotherapy (CT/RT). A significant decrease in risk to HNSCC was observed in the cases with deletion (CYP2A6*4B and CYP2A6*4C) or reduced activity genotypes (CYP2A6*9) of CYP2A6. This risk to HNSCC was further reduced significantly in tobacco users among the cases when compared to nontobacco users among the cases. The risk was also reduced to a slightly greater extent in alcohol users among the cases when compared to nonalcohol users among the cases. In contrast with decreased risk to HNSCC, almost half of the cases with variant genotypes of CYP2A6 (CYP2A6*1A/*4C+*1B/*4C+*4C/*4C and *9/*9) did not respond to the treatment. Likewise, the survival rate in cases receiving the treatment, after 55 months of follow-up was significantly lower in cases with deletion (6.3%) or reduced activity (11.9%) allele than in the cases with common alleles (41%). The present study has shown that CYP2A6 polymorphism significantly reduces the risk to HNSCC. Our data further suggested that CYP2A6 polymorphism may worsen the treatment outcome in the cases receiving CT/RT.

Deciphering Structural Alterations Associated with Activity Reductions of Genetic Polymorphisms in Cytochrome P450 2A6 Using Molecular Dynamics Simulations.

Cytochrome P450 (CYP) 2A6 is a monooxygenase involved in the metabolism of various endogenous and exogenous chemicals, such as nicotine and therapeutic drugs. The genetic polymorphisms in CYP2A6 are a cause of individual variation in smoking behavior and drug toxicities. The enzymatic activities of the allelic variants of CYP2A6 were analyzed in previous studies. However, the three-dimensional structures of the mutants were not investigated, and the mechanisms underlying activity reduction remain unknown. In this study, to investigate the structural changes involved in the reduction in enzymatic activities, we performed molecular dynamics simulations for ten allelic mutants of CYP2A6. For the calculated wild type structure, no significant structural changes were observed in comparison with the experimental structure. On the other hand, the mutations affected the interaction with heme, substrates, and the redox partner. In CYP2A6.44, a structural change in the substrate access channel was also observed. Those structural effects could explain the alteration of enzymatic activity caused by the mutations. The results of simulations provide useful information regarding the relationship between genotype and phenotype.

Genetic Variants in Smoking-Related Genes in Two Smoking Cessation Programs: A Cross-Sectional Study.

Previous studies have identified variants in genes encoding proteins associated with the degree of addiction, smoking onset, and cessation. We aimed to describe thirty-one single nucleotide polymorphisms (SNPs) in seven candidate genomic regions spanning six genes associated with tobacco-smoking in a cross-sectional study from two different interventions for quitting smoking: (1) thirty-eight smokers were recruited via multimedia to participate in e-Decídete! program (e-Dec) and (2) ninety-four attended an institutional smoking cessation program on-site. SNPs genotyping was done by real-time PCR using TaqMan probes. The analysis of alleles and genotypes was carried out using the EpiInfo v7. on-site subjects had more years smoking and tobacco index than e-Dec smokers (p < 0.05, both); in CYP2A6 we found differences in the rs28399433 (p < 0.01), the e-Dec group had a higher frequency of TT genotype (0.78 vs. 0.35), and TG genotype frequency was higher in the on-site group (0.63 vs. 0.18), same as GG genotype (0.03 vs. 0.02). Moreover, three SNPs in NRXN1, two in CHRNA3, and two in CHRNA5 had differences in genotype frequencies (p < 0.01). Cigarettes per day were different (p < 0.05) in the metabolizer classification by CYP2A6 alleles. In conclusion, subjects attending a mobile smoking cessation intervention smoked fewer cigarettes per day, by fewer years, and by fewer cumulative pack-years. There were differences in the genotype frequencies of SNPs in genes related to nicotine metabolism and nicotine dependence. Slow metabolizers smoked more cigarettes per day than intermediate and normal metabolizers.

Independent and Interactive Effect of CYPs and GSTs Genetic Variants and Tobacco Smoking on the Risk of Non-Small Cell Lung Carcinoma.

BACKGROUND: CYP and GST gene families detoxify tobacco carcinogens and have been linked to the risk of non-small cell lung carcinoma (NSCLC). AIM: Independent and combined effects of CYP and GST genetic variations and smoking on the risk of non-small cell lung carcinoma (NSCLC) and its sub-histological types. METHODS: We modelled an epistatic interaction via the effects of particular genotypes in two genes as OR (odds ratio), OR1, and OR2, a combination of both genotypes were characterized as ORcombine. In contrast, the two ORs' epistatic interaction for the individual genotypes has been represented as ORinteraction = ORcombine/(OR1 × OR2). RESULTS: The variant genotypes of CYP2A6 (OR:4.2, p <0.001), GSTT1 (OR:3.9, p <0.001), and GSTM1 (OR: 4.5, p <0.001) were showed a significant risk with NSCLC. GSTM1 (del.)/CYP2A6 (variant) genotype was associated with a higher risk of NSCLC (OR:12.5, p <0.001). GSTM1 (del.)/CYP2A6 (Ser/Pro+Pro/Pro) and GSTM1 (del.)/CYP2A13 (CT+TT) interacted redundantly (ORintraction = 0.66 and 0.64). A co-suppressive interaction was observed between GSTT1 (del.)/CYP2A6 (Ser/Pro+Pro/Pro) (ORintraction = 0.41). Simultaneously, both GSTT1/GSTM1 del. genotype was associated with a significantly higher risk to NSCLC. In contrast, GSTT1 del./GSTM1 del. genotype interaction displayed a co-suppressive effect (ORintraction = 0.15). CYP1A1(TC+CC)/CYP2A13(CT+TT)mutually interacted synergistically (ORintraction = 1.27).CYP1A1 (TC+CC)/GSTP1 (Val/Val+Ile/Val) genotype demonstrated an additive (ORintraction = 1) effect. GSTP1(Val/Val+Ile/Val) interacts with GSTT1 (del.) genotype exerted a suppressive effect (ORintraction = 0.69). CYP2A6 in smokers increased risk by 4.2 (p = 0.001) to 5.6 fold (p <0.001), while GSTM1 and GSTT1 were independent of smoking. CONCLUSION: Epistatic interactions revealed that CYPs/GSTs might follow a web of the interactions to modify the risk of NSCLC.

Genome-wide association study of letrozole plasma concentrations identifies non-exonic variants that may affect CYP2A6 metabolic activity.

OBJECTIVES: Letrozole is a nonsteroidal aromatase inhibitor used to treat hormone-receptor-positive breast cancer. Variability in letrozole efficacy and toxicity may be partially attributable to variable systemic drug exposure, which may be influenced by germline variants in the enzymes responsible for letrozole metabolism, including cytochrome P450 2A6 (CYP2A6). The objective of this genome-wide association study (GWAS) was to identify polymorphisms associated with steady-state letrozole concentrations. METHODS: The Exemestane and Letrozole Pharmacogenetics (ELPh) Study randomized postmenopausal patients with hormone-receptor-positive nonmetastatic breast cancer to letrozole or exemestane treatment. Germline DNA was collected pretreatment and blood samples were collected after 1 or 3 months of treatment to measure steady-state letrozole (and exemestane) plasma concentrations via HPLC/MS. Genome-wide genotyping was conducted on the Infinium Global Screening Array (>650 000 variants) followed by imputation. The association of each germline variant with age- and BMI-adjusted letrozole concentrations was tested in self-reported white patients via linear regression assuming an additive genetic model. RESULTS: There were 228 patients who met the study-specific inclusion criteria and had both DNA and letrozole concentration data for this GWAS. The association for one genotyped polymorphism (rs7937) with letrozole concentration surpassed genome-wide significance (P = 5.26 × 10-10), explaining 13% of the variability in untransformed steady-state letrozole concentrations. Imputation around rs7937 and in silico analyses identified rs56113850, a variant in the CYP2A6 intron that may affect CYP2A6 expression and activity. rs7937 was associated with age- and BMI-adjusted letrozole levels even after adjusting for genotype-predicted CYP2A6 metabolic phenotype (P = 3.86 × 10-10). CONCLUSION: Our GWAS findings confirm that steady-state letrozole plasma concentrations are partially determined by germline polymorphisms that affect CYP2A6 activity, including variants near rs7937 such as the intronic rs56113850 variant. Further research is needed to confirm whether rs56113850 directly affects CYP2A6 activity and to integrate nonexonic variants into CYP2A6 phenotypic activity prediction systems.